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1 .  INTRODUCTION 


The  purpose  of  this  study  is  to  investigate  new  techniques  for 
monitoring  density  fluctuations  in  a  turbulent  flow  field;  in  particular, 
a  new  three-level  gain  measurement  scheme  is  examined  with  the  aim  of  pro 
viding  important  information  about  the  turbulent  flow  field  behind  a  tur¬ 
ret  mounted  on  an  airplane.  For  this  application  there  are  certain  con¬ 
straints  under  which  the  density  diagnostic  is  to  function.  These  con¬ 
straints  include  a  spatial  resolution  of  approximately  1  cm,  and  the  capa 
bility  of  monitoring  1%  density  fluctuations  non-intrusively ,  instanta¬ 
neously,  accurately  and  remotely. 

Several  schemes,  which  may  be  used  to  monitor  species  concentra¬ 
tions  in  a  flow  field,  are  available,1  11  but  all  these  schemes  suffer 
from  one  or  more  drawbacks.  Methods  involving  shadowgraph,  schlieren, 
interferometry  and  resonance  absorption  of  light  can  only  provide  infor¬ 
mation  about  the  average  density  along  a  given  path  length.  Scattering 
techniques  based  upon  Raman  and  Rayleigh  processes,  in  principle,  can 
provide  the  necessary  spatial  and  temporal  resolution  needed  to  make 
such  measurements  in  a  turbulent  environment;  however,  Raman  scattering 
cross-sections  are  small  with  the  result  that  a  reasonably  good  signal 
to  noise  ratio  is  difficult  to  obtain,  while  the  Rayleigh  scattered  sig¬ 
nal  may  contain  significant  contributions  from  Mie  scattering.  Further¬ 
more,  under  combustion  conditions,  Robben11  has  concluded  that  for  the 
measurement  of  10%  root  mean  square  turbulent  fluctuations  or  less,  the 
Raman  scattering  technique  is  incapable  of  making  a  satisfactory  measure¬ 
ment. 

Recently,  Daily5  has  used  laser  induced  fluorescence  to  measure 
density  in  turbulent  reacting  flows.  An  important  feature  of  his  tech¬ 
nique  is  that  if  a  laser  is  used  to  saturate  a  particular  transition, 
the  fluorescence  signal  from  the  excited  state  is  independent  of  the 
quenching  rates.  The  possibility  of  using  this  technique  to  measure 
density  fluctuations  has  also  been  explored  in  this  report  in  Section  3. 
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The  main  emphasis  in  this  report  is  the  theoretical  evaluation 
of  a  three-level  gain  measurement  technique  as  a  possible  method  for  ob¬ 
taining  quantitative  information  about  density  fluctuations  in  an  air 
turbulent  boundary  layer.  Several  naturally  occurring  species  in  the 
atmosphere  can  be  used  to  make  such  a  measurement  and  details  about  three- 
level  gain  measurements  are  provided  in  Section  2.  In  Section  3,  tech¬ 
niques  for  monitoring  other  parameters  which  can  be  related  to  gas  den¬ 
sity  are  considered.  These  techniques  include  monitoring  laser  induced 
fluorescence  and  the  use  of  gas  breakdown  as  a  diagnostic  tool.  In  Sec¬ 
tion  4,  the  effect  of  turbulence  on  diagnostic  design  is  discussed,  and 
the  possibility  of  conducting  a  laboratory  experiment  to  test  the  diag¬ 
nostic  scheme  is  explored  in  Section  5.  Finally,  the  conclusions  are 
presented  in  Section  6. 
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2.  THREE-LEVEL  GAIN  MEASUREMENT  TECHNIQUES 


The  basic  three-level  gain  measurement  scheme  entails  using  a 
high  power  pulsed  laser  to  saturate  one  transition  of  a  molecule  and 
measuring  gain  on  a  coupled  transition  of  the  same  molecule  with  a  probe 
laser.  As  shown  in  Fig.  1,  the  pump  laser  saturates  the  transition  1  «-+■  2 
while  the  probe  laser  measures  gain  on  transition  2  «-*-  3.  The  power  of  the 
pump  laser  is  chosen  such  that  in  each  shot  the  transition  is  saturated 
even  though  there  may  be  variations  in  the  output  of  the  pump  laser  from 
shot  to  shot.  This  condition  is  realized  when  the  pump  laser  intensity 
is  always  greater  than  the  saturation  intensity  for  the  transition;  con¬ 
sequently,  pump  laser  variations  play  no  role. 

Saturating  transition  1  ■*-*■  2  means  that  a  specified  number  of 
molecules  (and  this  number  can  be  calculated)  will  be  excited  from  level 
1  to  level  2.  Measurement  of  small  signal  gain  on  transition  2  ■*-*■  3  can 
provide  quantitative  information  about  the  number  of  particles  in  level 
2.  This  in  turn  is  related  to  the  total  concentration  of  the  species 
through  the  saturation  condition.  Now,  if  the  concentration  of  the 
species  involved  in  a  gain  measurement  in  the  atmosphere  is  known,  we 
can  determine  the  atmospheric  density.  To  obtain  good  spatial  resolution 
required  for  making  measurements  in  a  turbulent  environment,  the  pump  and 
probe  lasers  are  placed  in  a  perpendicular  configuration  as  illustrated 
in  Fig.  2. 

A  list  of  naturally  present  species,  along  with  their  atmospheric 
concentrations,  is  given  in  Table  1.  It  can  be  seen  that  there  are  several 
species  which  might  be  suitable  for  a  gain  measurement.  Selection  of  a 
species,  however,  depends  upon  a  number  of  criteria  including  (a)  the 
availability  of  a  high  power  laser  having  a  lasing  frequency  resonant 
with  the  transition  frequency  of  interest  in  the  chosen  molecule,  (b) 
measurable  amount  of  gain  obtainable  on  the  coupled  transition,  and,  (c) 
availability  of  a  probe  laser  of  sufficient  intensity  to  measure  the 
gain.  There  are  very  few  high  power  lasers  available  in  the  visible  and 
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Fig.  2  Arrangement  of  Pump  and  Probe  lasers  in  the  three  level 
scheme  to  make  spatially  resolved  measurements. 


TABLE  1. 


Concentration  of  Gases  in  Air 


Constituent 

ppm  by  Volume 

N2 

7.  808  x  10^ 

°2 

2.  09  x  105 

5jC 

h2° 

2.  5  x  10^ 

A 

9.  34  x  103 

C°2 

330 

Ne 

18.  18 

He 

5.  24 

CH4 

1.  6 

Kr 

1.  14 

H2 

0.  5 

n2o 

0.  28 

Xe 

0.  087 

CO 

0.  075 

For  tropical  atmospheres  at  sea  level, 
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UV.  This  means  that  saturating  an  electronic  transition  is  perhaps  not 
practical.  Ideally,  one  would  like  to  work  with  nitrogen  or  oxygen  be¬ 
cause  these  gases  have  appreciable  concentrations  in  the  atmosphere. 
However,  there  are  no  laser  frequencies  which  are  resonant  with  the  two 
coupled  allowed  transitions  in  these  molecules. 

The  amount  of  wate~  vapor  in  the  atmosphere  varies ,  and  for  tro¬ 
pical  atmospheres  at  sea  levex,  the  water  vapor  concentration  is  2.5% 
by  volume.  The  water  molecule  is  one  of  the  molecules  which  has  been 
investigated  in  this  study  for  a  possible  three-level  gain  measurement. 
Another  molecule  considered  is  C02 ,  which  has  a  concentration  of  .03%  by 
volume.  Argon  and  other  inert  gases  have  not  been  studied  because  elec¬ 
tronic  transitions  are  involved.  Methane,  N^O  and  other  trace  gases  are 
possible  candidates  for  a  three-level  gain  measurement,  but  their  con¬ 
centration  in  the  atmosphere  is  too  small  to  produce  an  appreciable  gain 
necessary  for  an  accurate  measurement  of  density. 

2 . 1  Three-Level  Gain  Measurement  in  CO^ 

In  the  lowest  electronic  state  of  C02  there  are  several  coupled 
transitions  involving  vibrational  levels  which  can  be  employed  for  a 
three-level  measurement.  However,  the  suitability  of  a  particular  tran¬ 
sition  for  making  a  density  measurement  depends  upon  the  availability  of 
lasers.  With  the  use  of  the  existing  commercially  available  laser  sys¬ 
tems,  the  three-level  gain  measurement  in  CO 2  may  be  made  in  the  follow¬ 
ing  three  ways:  (i)  2.7  ym  HF  laser  pumping  of  the  000  +-*■  021  or  000  «-*011 
transition  and  gain  measurement  on  the  021  *-*■  020  transition  with  a  4.3  ym 
CW  diode  laser;  (ii)  4.3  ym  HBr  laser  pumping  of  the  000  «-*■  001  transition 
and  gain  measurement  on -the  001  *-*■  100  transition  with  a  10.6  ym  CO  probe 
laser;  and,  (iii)  pumping  of  the  000  •*->■  001  transition  with  a  4.3  ym 
HBr  laser  and  gain  measurement  on  the  001  000  transition  with  a  4.3  ym 

CW  diode  laser. 

Detailed  theoretical  calculations  have  been  made  and  these  calcu¬ 
lations  show  that  scheme  (i)  is  the  most  promising  of  the  three  schemes 
outlined  above.  Consequently,  in  the  following  we  will  discuss  scheme  (i) 
in  detail,  while  the  other  two  schemes  will  be  described  only  briefly. 
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The  vibrational  levels  021  and  Oil  of  CO^  are  in  Fermi  resonance. 
Therefore,  we  assume  that  these  two  levels  are  in  equilibrium  at  all  times. 
Several  HF  laser  lines  are  in  near  resonance  with  different  rotational 
lines  of  combination  bands  000  021  and  000  +-»•  Oil  in  C02-  The  selection 

of  a  particular  HF  laser  line  is  based  upon  the  following  criterion: 

(a)  The  HF  laser  line  should  be  in  resonance  with  a  high 
J  rotational  line  on  an  R-branch  transition.  Theore¬ 
tical  calculations  given  in  Appendix  II  show  that 
when  a  high  intensity  laser  is  resonant  with  a  high 

J  rotational  line  in  an  R-branch  transition,  it  is 

possible  to  excite  more  than  half  of  the  ground  state 

* 

particles  to  the  excited  vibrational  manifold; 

(b)  The  HF  laser  line  selected  for  pumping  must  be  one 
of  the  stronger  lines  of  the  HF  laser  to  achieve 
the  necessary  saturation  condition; 

(c)  The  laser  intensity  needed  to  cause  saturation  should 
be  such  that  it  must  not  cause  breakdown  in  the  air. 

After  taking  into  consideration  the  above  criterion,  it  turns  out 
that  the  P(4)  line  of  the  transition  2  •*-*•  1  in  an  HF  laser  is  best  suited 

to  saturate  the  000  +-*■  021  transition  in  CO..  The  P(4)  line  (V  =  3622.663 

-1  2  P-1 
cm  )  is  in  near  resonance  with  the  R(12)  line  (V  =  3622.492  cm  )  of  the 

P 

100  ■<-*■  021  transition  in  C02.  The  exact  frequencies  of  HF  laser  lines  are 
taken  from  Ref.  12  while  the  frequency  and  the  strength  of  rotational 
lines  in  a  particular  band  of  C02  are  obtained  from  a  computer  program 
which  uses  the  parameters  compiled  in  Ref.  13. 


* 

To  be  published  in  Optics  Letters,  December  1979. 


The  laser  intensity  required  to  saturate  a  transition  can  be  cal¬ 
culated  from  Eq.  (AI-10)  derived  in  Appendix  I.  To  calculate  saturation 
intensity  one  needs  the  characteristic  time  for  vibrational  relaxation 
(T  ) ,  the  Einstein  coefficient,  (A)  ,  for  vibration-rotation  transitions 
corresponding  to  the  R(12)  line  of  the  000  *-*■  021  transition  in  CO^  and 
the  width  (Av  )  of  the  rotational  line  R(12) .  The  vibrational  relaxation 

Li 

time  constant  (T  )  for  level  021  has  been  taken  to  be  the  characteristic 
v 

time  for  the  V-V  transfer  process  which  is  the  fastest  rate  of  deactiva¬ 
tion  of  the  021  level  in  CO^. 

C02  (021)  +  N2  (v  =  0)  -*■  C02  (020)  +  N  (v  »  1)  . 

The  rate  constant  for  this  reaction  is14  1.6  x  106  sec  \  Consequently, 

=  6  x  10  sec.  Other  deactivation  rates  for  the  excited  level  021 

are  much  smaller  and  the  value  6  x  10  7  sec  should  be  correct  to  within 

a  factor  of  2  of  the  true  value  of  T  . 

v 

The  Einstein  coefficient  for  the  particular  vibration-rotation 
transition  has  been  calculated  from  the  following  relation: 


A 


(2.48  x  lO1^)  exp  (V 


1.4387/T)  , 


(1) 


where 


S  =  Line  strength  in 


cm 


,  cm 


-2 


molecule  obtained  from  a 


computer  program  based  upon  parameters  of  Ref.  13. 


A  =  Transition  wavelength  in  )im. 


-1 


V  =  Energy  of  the  lower  level  in  cm 

1j 


For  the  R(12)  line  of  000  021  transition  in  CO^, 


T  =  300  K,  A  =  2.76  pm. 

Substituting  these  values  into  Eq.  (1)  we  find  A  =  9.38  sec  \ 

The  width,  Av  ,  of  the  rotational  line  R(12)  has  been  taken  to 

Xj 

be  the  mean  full  width  at  half  maximum  of  the  CO  molecule  at  1  atmo- 

sphere  pressure  as  listed  in  Ref.  13  (Av  =  0.07  cm  ) . 

L 

Substituting  the  value  of  x  ,  A,  and  Av  in  Eq.  (AI-10) ,  we  find 

V  L 

that  the  laser  intensity  required  to  saturate  the  R(12)  line  in  CO  must 

,  2  2 
be  much  greater  than  0.185  MW/cm  .  To  saturate  the  vibrational  transi¬ 
tion  000  ■*->■  021  completely,  however,  this  intensity  must  be  divided  by 
the  Boltzmann  factor  for  the  R(12)  line  (see  Appendix  I).  Taking  into 
account  the  Fermi  resonance  between  levels  021  and  Oil,  the  Boltzmann 
factor  for  the  R(12)  line  is  0.0173.  Therefore,  the  laser  intensity 

needed  to  saturate  the  vibrational  band  000  •<-+•  021  must  be  much  greater 
2  2 

than  11  MW/cm  .  Thus,  tens  of  MW/ cm  of  laser  power  will  be  required 
to  saturate  the  combination  band  000  ■*-*■  021. 

A  laser  intensity  of  tens  of  megawatts  per  square  cm  at  2.7  pm 
does  not  cause  breakdown  in  clean  air.  Even  when  the  air  is  mixed  with 
particulate  matter,  there  will  be  no  breakdown  in  air  at  2.7  pm  if  the 
laser  pulse  duration  is  1  psec  or  shorter.  The  choice  of  the  laser 
pulse  duration,  t  ,  depends  upon  (T  _<  Tv) .  Since  Ty  =  .6  psec,  the 
whole  experiment  will  be  completed  in  less  than  1  psec.  Thus,  under 
experimental  conditions  of  interest  to  us,  air  breakdown  should  not  oc¬ 
cur. 
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S  =  1.035  x  10 
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2 . 2  Calculation  of  the  Absorption  Coefficient 

The  knowledge  of  the  absorption  coefficient  in  air  at  the  pump 
laser  frequency  is  necessary  to  determine  how  far  along  the  path  the 
laser  beam  will  be  able  to  saturate  the  gas.  The  absorption  coefficient 
in  cm  3  has  been  calculated  using  the  following  relation: 


K(V) 


S  Av  N 
L 


2  IT 


(2 


where 

S  =  Line  strength  in  (cm  molecule  . 

N  =  Number  density  of  absorbers,  cm  3. 

Vp  =  Pump  laser  frequency,  cm  3 . 

=  Line  center  frequency  of  the  transition,  cm  ^ . 

Av  =  Full  width  at  half  maximum,  cm  \ 

L 

—18  —1 

For  the  HF  laser  line  P(4)  2  -*■  1,  S  =  1.035  x  10  cm  molecule  , 

N  =  7.35  x  lO^  molecules/cm3  at  sea  level,  (V  -  V  )  =  0.171  cm  3 

-1  P  ° 

and  Av  =  0.07  cm  .  Substituting  these  values  in  Eq.  (2),  we  get 
L  -4  -1 

<(v  )  =  1.71  x  10  cm  .  Thus,  the  HF  laser  will  be  able  to  satu- 
P 

rate  the  000  021  combination  band  in  CO^  over  a  path  length  of 

several  meters  in  air. 

2.3  Calculation  of  Small  Signal  Gain 


Thus  far  we  have  concentrated  on  the  calculation  of  pump 
laser  parameters  such  as  laser  intensity,  pulse  duration  and  laser 


absorption  coefficient  in  air  required  to  saturate  the  000  •*-»-  021  com¬ 
bination  band  in  CO^.  By  selecting  an  appropriate  rotational  line  we 
have  also  tried  to  maximize  the  number  of  particles  excited  to  the  021 
level  in  CO Now  we  concentrate  on  the  calculation  of  probe  laser 
parameters  and  on  the  maximization  of  gain  on  the  coupled  transition 
021  +-+  020. 

In  a  thermalized  rotational  distribution,  the  gain  on  a  P-branch 
transition  is  larger  than  on  Q-  or  R-branch  transitions.  Consequently, 
for  a  given  rotational  partition  function,  the  maximum  gain  is  obtained 
when  the  rotational  level  with  the  most  particles  in  the  excited  state 
is  involved  in  the  P-branch  transition.  On  the  coupled  transition 
021  ■*-*■  020  in  CC^,  the  P  (16)  line  with  frequency  2314.22  cm  3  will  give 
the  maximum  gain.  To  measure  this  gain,  a  CW  diode  laser  will  be  used. 

This  laser  can  be  fine  tuned  to  match  exactly  the  P(16)  line  of  C02  on 
the  021  «-*  020  transition. 

The  small  signal  gain,  G,  on  the  P(16)  line  of  the  021  •*-*  020 
transition  in  C02  has  been  calculated  using  the  following  expression: 

- _  (3) 

4TT2  AVl  1  +  4  (■  La-  °) 

where  An  is  the  population  inversion  on  the  P(16)  line,  A  is  the  Ein- 

J  tJ 

stein  coefficient  for  spontaneous  emission,  V  is  the  frequency  of  the 

L 

probe  laser,  and  the  other  symbols  have  been  defined  previously. 

The  population  inversion  AN  is  calculated  assuming  that  satu- 

J 

ration  with  the  pump  laser  will  place  approximately  one  third  of  the 

total  molecules  in  the  021  level  because  there  is  Fermi  resonance  between 

levels  021  and  011.  Substituting  the  values  An^  =  7.35  x  1015  x  .33  x  .0369 

=  8.95  x  103-3  molecules/cm3,  A,  =  173  sec  3,  Av_  =  0.07  cm  3,  (VT  -  V  )  s  0 

u  L  ^  L  o 

and  X  =  4.32  ym  in  Eq.  (3) ,  we  get  G(V  =  4.32  ym)  =  0.017  cm 

L 


G  (VL)  = 


-12- 


Measurement  of  gain  at  2314.22  cm  1  (4.32  pm)  may  become  difficult 

if  there  is  appreciable  absorption  of  the  probe  laser  in  air.  The  4.3  ym 

band  in  CO  originating  from  the  ground  state  is  very  strong.  However,  the 
z  -1 

frequency  2314.22  cm  is  appreciably  shifted  from  any  of  the  vibration- 
rotation  lines  in  the  CO^  000  ■*-»■  001  band.  To  see  how  much  power  the  probe 
laser  will  lose  as  a  consequence  of  this  absorption,  we  calculate  the  ab¬ 
sorption  coefficient  in  air  at  2314.22  cm  1  from  Eq.  (2) .  This  absorption 

-4  -1 

coefficient  comes  out  to  be  7.5  x  10  cm  .  Thus,  the  probe  laser  can 
traverse  a  distance  of  several  meters  without  appreciable  attenuation  in 
the  air. 

In  summary,  therefore,  the  three-level  gain  measurement  technique 
involving  the  021  ■«-»■  000  combination  band  in  CO^  yields  a  gain  of  0.017  cm 
with  the  use  of  commercially  available  laser  systems.  In  Section  IV,  we 
will  see  whether  a  gain  of  0.017  cm  ^  is  sufficient  to  provide  quantita¬ 
tive  information  about  turbulent  density  fluctuations  in  atmosphere.  For 
convenience,  first  we  summarize  in  Fig.  3  the  salient  features  of  the  4.3 
ym  gain  measurement  on  021  -*-*■  000  transition. 

Now  let  us  briefly  consider  other  techniques  for  making  a  three- 

level  gain  measurement  in  CO^.  One  scheme  involves  pumping  the  000  ■*-*■  001 

transition  in  CO^  with  a  4.3  ym  HBr  laser  and  a  gain  measurement  on  the 

OOl^lOO  transition  with  a  10.6  ym  CC>2  probe  laser.  The  pump  and  probe 

laser  requirements  and  other  salient  features  of  this  scheme  are  summarized 

in  Fig.  4.  Calculations  of  the  various  parameters  have  been  made  following 

the  procedure  outlined  above.  From  Fig.  4,  it  can  be  seen  that  one  can 

-4  -1 

only  get  a  gain  of  2  x  10  cm  which  is  much  smaller  than  the  gain  of 
0.017  cm  1  obtained  from  021  «->•  020  transition  considered  above. 

Another  three-level  gain  measurement  in  CO^  involves  pumping  the 
000  ■<-*■  001  transition  with  a  4.3  ym  HBr  laser  and  a  gain  measurement  on 
the  same  001  «-*•  000  transition  with  a  4.3  ym  CW  diode  laser.  This  tech¬ 
nique  is  based  upon  the  fact  that  under  certain  conditions  (which  are 
specified  in  Appendix  II) ,  optically  pumping  an  R-branch  transition  at 
high  J  number  leads  to  gain  at  those  frequencies  which  are  smaller  than 
the  pump  laser  frequency.  For  this  scheme,  the  pump  and  probe  laser  re¬ 
quirements  are  given  In  Fig.  5.  It  can  be  seen  that  this  scheme  is  not 


-13- 


Fermi 

—  Vibration  Level 

02  1 

Resonance 

=  10  1 

— 

- s-  Jb15 

J  -  13  = 

?=--■■■  X  P  (16) 

R(12)  / 

X^o- 

4.32  /am 

VQ  ■  2.76  /um  / 

J  -  16 

Type 

Laser  frequency 

Absorption/Gain  line 

Saturation  intensity 

Pulse  duration 

Laser  fluence 

Atmospheric  absorption 
coefficient 

Gain  at  4.  32  pm 


Pump  Laser 
Characteristics 

Probe  Laser 
Characteristics 

HF  Pulsed  Laser 

CW  Diode  Laser 

2.  76  pm,  (2  h>  1,  P(4) ) 

4.  32  pm 

CC>2,  000  «-*  021,  R(  1 2) 

C02,  021  020,  P(16) 

2 

>  ]  0  MW  /  cm 

— 

s  0.  6  psec 

CW 

2 

>  6  J/cm 

— 

,  -4  -1 

.4 

1.  71  x  10  cm 

7.  5  x  10 

- 1 

0.  0 17  cm 

Fig.  3  Salient  features  of  4.  3  pm  gain  measurement  on  hot  band 
(021  «-»  020)  in  CO,. 
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suitable  for  our  application  due  to  two  reasons.  First,  the  gain  obtain¬ 
able  on  the  001  -*-►  000  transition  is  only  7.67  x  10  3cm_1,  which  is  smaller 
than  the  gain  of  0.017  cm  1  obtainable  from  the  021  ■*-+  020  transition  con¬ 
sidered  above.  Secondly,  the  probe  laser  absorption  coefficient  in  the 
air  is  very  high  (  =  0.46  cm  ) ,  with  the  result  that  most  of  the  energy  in 
the  probe  laser  will  be  absorbed  in  the  atmosphere  by  the  time  the  probe 
laser  reaches  the  detector. 

2 . 4  Three-Level  Gain  Measurement  in  Water  Vapor 

The  amount  of  water  vapor  in  the  atmosphere  varies  with  the  al¬ 
titude,  geographic  location  and  the  time  of  the  year.  For  tropical  at¬ 
mospheres  at  sea  level,  the  water  vapor  concentration  is  2.5%  by  volume.13 
This  concentration  is  approximately  8  times  the  concentration  of  CO ^  in  the 
atmosphere  at  sea  level  and  justifies  studying  the  water  molecule  to  explore 
the  poss.ibility  of  a  three-level  gain  measurement  in  H^O  with  higher  gain. 

For  a  possible  three-level  gain  measurement  in  H^O,  two  coupled 
transitions  have  been  identified.  These  are  the  4^  •<-*■  4  transition  on 
the  001  ►  000  band,  and  the  6^  5^  transition  on  the  001  -*-*•  020  band. 

The  HF  laser  line  1  0,  P(5),  is  in  near  resonance  with  the  4  4 

transition  of  the  001  000  band  in  H^O.  The  gain  may  be  measured  on  the 

6^3  *-*■  5j.q  transition  on  the  001  -*-*  020  band  either  with  a  CW  diode  laser 
or  with  a  CW  H^O  laser.  The  salient  features  of  this  scheme  are  presented 
in  Fig.  6.  Here  several  points  should  be  noted.  First,  the  calculated 
gain  at  28  ym  turns  out  to  be  slightly  larger  than  the  4.3  ym  gain  measu¬ 
rement  in  CO^  (Fig.  3) .  However,  this  gain  is  only  for  tropical  atmos¬ 
pheres  at  sea  level.  As  one  moves  to  other  model  atmospheres,  the  gain 
decreases.  Furthermore,  at  higher  altitudes,  the  concentration  of  H^O 
in  air  decreases  and  so  does  the  gain.  For  instance,  as  shown  in  Fig.  7, 
at  an  altitude  of  6  Km,  we  can  only  obtain  a  gain  of  0.001  cm  1  for  tro¬ 
pical  atmospheres.  Thus,  this  scheme  becomes  increasingly  unattractive 
as  the  altitude  increases. 

The  other  drawbacks  of  this  scheme  are:  (i)  the  gain  is  measured 
at  longer  wavelengths  (28  ym)  where  detector  sensitivity  is  lower  yet;  and 
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Fig.  6  Salient  features  of  28  pm  gain  measurement  in  H  O. 
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altitude  for  five  different  model  atmospheres 


(ii)  the  laser  intensity  required  to  saturate  the  431  trans:*-t;’-on  in 

the  001  -<->■  000  band  is  rather  large  and  this  intensity  may  cause  breakdown 
in  the  air. 

A  summary  of  results  for  the  three-level  gain  measurement  schemes 
are  presented  in  Table  2  for  comparison.  It  is  clear  from  this  table  that 
the  most  promising  of  all  the  three-level  schemes  is  the  one  which  invol¬ 
ves  pumping  the  R(12)  line  of  the  021  +-*■  000  band  in  CO^  with  an  HF  laser 
and  the  gain  measurement  on  the  P(16)  line  of  the  021  020  band  with  a 

CW  diode  laser.  In  the  next  section,  we  will  discuss  those  schemes  which 
monitor  parameters  other  than  gain  to  yield  information  about  turbulent 
density  fluctuations. 


-20- 


warn 


TABLE  2 


3.  OTHER  DENSITY  FLUCTUATION  DIAGNOSTIC  TECHNIQUES 


In  this  section,  we  consider  those  techniques  which  monitor  para¬ 
meters  other  than  small  signal  gain  to  obtain  information  about  the  den¬ 
sity  fluctuations  in  a  turbulent  boundary  layer.  In  this  category,  two 
schemes  have  been  examined.  The  first  is  the  use  of  laser  induced  air 
breakdown  as  a  diagnostic  tool  and  the  second  is  the  ultra  violet-visible 
fluorescence  probes  to  monitor  turbulent  density  fluctuations.  Both  these 
schemes  are  discussed  below. 

3 . 1  The  Use  of  Laser  Induced  Air  Breakdown  as  a  Diagnostic  Tool 

The  breakdown  of  a  gas  from  the  radiation  of  a  focused  TEA  laser 

is  now  a  common  laboratory  phenomenon,  occurring  at  a  range  Of  different 

1 5  16 

laser  wavelengths,  pressures,  intensities  and  gases."  '  Many  of  the 
studies  of  this  breakdown  have  been  concerned  with  the  complex  mechanisms 
that  initiate  the  breakdown  by  supplying  the  initial  electrons  that  are 
accelerated  by  the  field  in  the  cascade  process.  This  breakdown  of  the 
gas  may  be  broken  up  into  several  temporal  regimes.  The  initiation  of 
the  breakdown  is  dependent  upon  the  presence  of  at  least  one  electron  with¬ 
in  the  focal  volume.  The  multiphoton  processes  responsible  for  the  cre¬ 
ation  of  this  1  lucky  electron"  are  still  a  subject  of  both  theoretical 
and  experimental  investigation,  and  the  differences  between  breakdown  in 
"clean"  and  "dirty,"  or  laboratory  air  are,  in  large  part,  explained  by 
the  presence  of  impurities  capable  of  easily  donating  electrons  to  start 
the  cascade  process.  Once  present,  an  electron  will  be  accelerated  by 
the  intense  field  in  the  focal  volume  of  the  laser  to  gain  sufficient 
energy  to  ionize  a  neutral  particle.  This  results  in  the  production  of 
two  low  energy  electrons  that  are  again  subsequently  accelerated  in  the 
field.  This  process  repeats  itself  leading  to  a  high  degree  of  ioniza¬ 
tion  in  an  extremely  short  time.  In  the  analysis  that  follows,  we  shall 
be  concerned  only  with  the  cascade  process  rather  than  multiphoton  pro¬ 
cesses  that  may  be  of  importance  during  the  early  stages  of  ionization. 
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Fig.  8  Summary  of  observations  for  a  pressure  of 
order  1  atm.  Vertical  scale  for  I  and  P  is 
arbitrary.  This  figure  has  been  taken  from 
Ref.  17. 


In  Fig.  8  we  see  the  qualitative  behavior  of  air  breakdown  at  a  pressure 

17  18 

of  nearly  one  atmosphere.  '  For  higher  laser  fluxes,  an  order  of  magni- 

19 

tude  or  more  above  threshold,  and  for  clean  N2,  the  breakdown  is  one  to 
two  orders  of  magnitude  more  rapid.  The  purpose  of  the  following  investiga 
tion  will  be  to  examine  the  dependence  of  electron  density  upon  neutral 
density,  time,  and  laser  flux  during  the  initial  stages  of  breakdown. 

In  Fig.  9  we  see  a  possible  arrangement  by  which  plasma  breakdown 

and  a  measurement  of  the  time  to  opacity  or  complete  ionization  can  give 

information  on  the  spectral  density  and  frequency  of  density  fluctuations. 

At  this  point  we  will  merely  note  the  well  known  correlation  between 

spatial  measurements  of  a  fluctuating  quantity  at  a  single  instant  of 

time  and  the  temporal  fluctuation  (the  actual  quantity  desired)  at  one 
19 

point  in  space.  In  our  proposed  experimental  configuration,  we  attempt 
to  remove  the  one  disadvantage  associated  with  using  laser  breakdown  to 
obtain  quantitative  information  on  local  properties,  and  that  is  the 
shot  to  shot  variation  of  the  TEA  laser  itself. 

In  the  configuration  indicated,  the  pulse  from  the  TEA  laser 
passes  through  a  50%  beamsplitter.  At  each  focal  point,  a  breakdown 
will  occur,  and  the  degree  of  ionization  will  be  such  as  to  severely  at¬ 
tenuate  the  CW  C02  probe  laser  passing  through  these  focal  points.  It  is 
the  "time  to  opacity"  measured  at  detectors  1  and  2  that  will  provide  in¬ 
formation  on  the  local  state  of  the  gas  in  the  breakdown  volume.  The 
measurements  must  be  made  as  we  shall  see  below,  in  the  order  of  tens  of 
nanoseconds,  which  will  still  ensure  that  we  are  concerned  with  small,  i.e. 
"point"  volumes  of  gas.  With  additional  great  attention  placed  upon  ob¬ 
taining  high  beam  quality,  this  will  minimize  the  differences  in  intensity 
between  focal  points  1  and  2.  Thus,  any  shot  to  shot  variation  in  the 
pulse  will  manifest  itself  in  the  same  manner  at  each  breakdown  region. 

A  stable  CW  probe  laser,  also  CC>2,  passes  through  the  other  side  of  the 
beamsplitter  as  shown.  This  will  cause  the  probe  laser  to  pass  through 
each  focal  region.  If  the  second  mirror  and  lens,  as  well  as  the  second 
detector  were  on  an  optical  rail,  then  the  second  focal  point  could  be 
moved  easily  to  measure  the  spatial  correlations  of  the  breakdown  times 
without  a  need  to  realign  either  laser.  Each  detector  will  monitor  the 
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Fig.  9  Experimental  configuration  for  density  measurement 
using  gas  breakdown. 
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time  required  to  reach  opacity,  and  it  will  be  shown  how  the  difference  in 
time  required  to  reach  full  ionization  at  focal  points  1  and  2  will  be  a 
sensitive  measure  of  the  difference  in  local  neutral  particle  density.  In 
order  to  insure  that  the  detectors  sample  only  the  CW  radiation,  we  can 
arrange  that  the  CW  laser  has  polarization  different  from  that  of  the  pulsed 
laser,  and  that  a  suitable  polarizer  be  placed  before  each  detector.  The 
following  order  of  magnitude  analysis  will  provide  an  estimate  of  the  time 
required  to  reach  opacity,  (i.e.,  the  density  associated  with  the  local 
plasma  frequency)  and  the  sensitivity  of  this  rise  time  to  the  local  par¬ 
ticle  density. 

3.1.1  Order  of  Magnitude  of  the  Time  Required  To  Reach  Opacity 


The  plasma  frequency  U)  depends  upon  the  number  of  free  electrons 
present  that  interact  with  radiation  of  a  given  frequency.  When  the  plasma 
frequency  is  above  or  of  the  order  of  the  laser  frequency,  significant  ab¬ 
sorption  results.  The  time  to  reach  this  state  is  what  we  have  called 

4  t 

above  the  "time  to  opacity."  Thus,  cj  .=  5.64  x  10  (n  )  ,  where  w  is 

3  P  ®  P 

in  rad/sec,  and  n  in  electrons/cm  .  For  10.6  micron  radiation,  U)=  2ttc/A  = 

14  6  19  3 

1.77  x  10  rad/sec.  For  U)  «  w,  n  a  10  electrons/cm  .  Thus,  at  one 

p  e 

atmosphere,  complete  ionization  produces  an  electron  density  that  cor¬ 
responds  to  the  plasma  frequency.  The  time  to  yield  this  required  number 
of  electrons,  and  the  dependence  of  the  electron  concentration  upon  tne 
initial  neutral  density  will  be  estimated  in  the  following. 

We  consider  first  the  full  energy  balance  describing  the  energy 
increase  of  the  electrons 
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where  is  the  free-free  absorption  coefficient  for  neutral  particles,  I 

is  the  laser  flux ,  and  K.  is  the  ion  absorption  coefficient.  The  electron 

1  * 

density  is  given  by  n^;  n^  is  the  excited  level  number  density  with  ener¬ 
gy  E  ;  n  is  the  neutral  density;  I  is  an  effective  ionization  potential; 
a  p 

m  is  the  mass  of  an  electron  and  M  is  the  mass  of  a  neutral.  We  neglect 
compared  with  K  ,  true  certainly  until  the  ionization  becomes  appreci¬ 
able.  Further, 
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where  f  is  the  neutral-electron  collision  frequency.  Considering  only 
processes  that  create  electrons  in  Eq.  (4),  ve  may  write,  as  an  order 
of  magnitude  approximation. 
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From  this  the  number  of  collisions  required  to  raise  the  energy  of  an 
electron  by  hv  is  approximately 


TTmc 

I 


(7) 


Since  we  wish  the  avalanche  process  to  be  relatively  slow,  we  choose  a 

9  2 

threshold  laser  intensity  of  I  =  2.8  x  10  watts/cm  .  For  a  10.6  micron 
photon,  V  =  2.8  x  1033  and  for  the  intensity  cited  above,  we  find  approxi¬ 
mately  2  collisions  are  needed  for  an  energy  increase  by  hv.  Thus,  if  we 
take  as  the  air  constituent  with  the  lowest  ionization  potential  (12.3 
and  we  neglect  NO  (ionization  potential  9.6  ev)  the  number  of  collisions 


ev)  , 
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required  to  raise  a  neutral  particle  to  the  ionization  potential  is  given 

by  Eq.  (7),  multiplied  by  12.3  ev/hv  100,  or  200  collisions  necessary 

22 

to  ionize  an  O  molecule.  Using  data  from  Phelps  and  Hake  for  O  ,  we 

2  10  o  2 

find  that  f  «  8.4  x  10  collisions/sec.  at  T  s  300  K.  Thus,  200  (col- 

10  60  -9 

lisions/electron  pair)/  8.4  x  10  (collisions/sec . )  ~  2.3  x  10  sec,  the 
time  to  create  one  electron  pair.  If  we  consider  the  avalanche  to  begin 
when  we  have  reached  109  -  10 9-0  electrons/cm3,  (see  Fig.  8)  then  we  must 
estimate  the  time  required  to  increase  this  original  number  of  electrons 

to  an  electron  density  corresponding  to  the  plasma  frequency  for  10.6 

19  3  9 

micron  radiation,  i.e.,  10  electrons/cm  .  If  we  begin  with  10  prim¬ 
ing  electrons/cm3  then  the  number  of  generations  of  electrons  is  given 
n  19  9 

as  2  a  10  /10  .  For  this  case  n  =  33  generations  of  electrons.  If  we 

12 

begin  with  10  electrons,  then  the  number  of  generations  needed  to  pro- 
19 

duce  10  electrons  is  23.  Thus,  between  23  and  33  generations  of  elec¬ 
trons  are  needed  to  reach  full  ionization  with  the  above  intensity  and 
for  10.6  micron  radiation.  With  one  generation  of  electron  pairs  re¬ 
quiring  approximately  2  ns,  then  we  are  considering  the  time  to  breakdown, 
again  for  the  conditions  specified  here,  of  the  order  of  45-65  nsec.  It 
should  be  noted  that  use  of  plasma  opacity  as  a  fast  "shutter"  requires 

an  intensity  two  orders  of  magnitude  larger  than  specified  here,  and 

18 

consequently,  the  rise  time  is  sub  nanosecond.  We  thus  estimate  that 
the  experiment  proposed  will  occur  on  a  time  sale  of  the  order  of  several 
tens  of  nanoseconds,  well  within  the  reach  of  a  modern  oscilloscope.  A 
more  quantitative  estimate  of  the  dependence  of  this  rise  time  on  the 
neutral  density  will  be  conducted  in  the  following  subsection. 


3.1.2  Dependence  of  Rise  Time  on  Neutral  Densitv 


Variations  in  the  output  power  of  the  laser  causing  the  break¬ 
down,  as  well  as  variations  in  the  neutral  density  of  the  medium  in  which 
the  breakdown  occurs,  will  cause  the  time  to  reach  opacity  to  vary.  We 
will  also  see  that  the  dependence  upon  initial  electron  density  (assum¬ 
ing  that  the  threshold  condition  is  satisfied)  will  be  relatively  slight. 

That  is,  as  noted  above,  a  change  in  the  number  of  priming  electrons/cm3 

12  7  ... 

from  10  to  10  only  changes  the  time  to  reach  complete  ionization  by 

a  factor  of  two. 
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We  now  consider  the  equation  for  electron  conservation,  in  stan¬ 


dard  notation, 


dn 
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Bn2  n.  . 
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Since  we  are  considering  processes  that  will  produce  the  fastest  possible 
rise  time,  we  neglect  all  loss  terms  in  the  energy  equation,  Eq.  (4).  We 
will  do  the  same  for  the  continuity  equation,  although  an  iterative  pro¬ 
cedure  on  the  loss  term  in  the  electron  continuity  equation  could  easily 
be  employed.  With  this  simplification,  the  energy  equation  for  electrons 
becomes 
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where  f  (T  )  is  a  function  of  time.  In  principle,  the  laser  flux,  I,  may 
also  be  a  function  of  time.  With  the  neglect  of  loss  processes  to  es¬ 
tablish  a  lower  bound  on  the  time  to  opacity,  certainly  valid  until  an 
appreciable  fraction  of  the  neutral  particles  are  ionized,  we  uncouple 

the  electron  energy  and  continuity  equations.  We  assume  that  the  effec- 

.  22 
tive  collision  frequency  is  proportional  to  a  mean  electron  velocity 

and  that  f  ~  (kT  )  .  For  convenience,  we  will  also  choose  a  reference 

e 

frequency  and  temperature  to  enable  us  to  write 
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Eq.  (7)  may  now  be  written  as 


This,  in  principle,  allows  consideration  of  pulse  shapes  upon  the  increase 
of  the  electron  temperature.  In  the  following,  however,  we  will  restrict 
the  analysis  to  the  case  in  which  the  pulse  is  time  independent.  There¬ 
fore  , 
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with  f  ~  n  ,  we  see  that  as  soon  as  the  electron  temperature  is  notably 
R  a 

above  the  initial  electron  temperature, 
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We  now  investigate  the  dependence  of  the  electron  density  upon  electron 
temperature.  Loss  terms  in  the  electron  density  equation  are  omitted 
and  we  note  that 


c. 
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Jj 

with  o  ~  T  and  <  v  >  ~  T  .  Again,  choosing  a  reference 

i  “  0  6  0 

Te^  (not  necessarily  the  same  as  previously) ,  we  write 


temperature. 
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The  equation  for  the  increase  of  electron  density  may  now  be  written  as 
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The  strong  non-linearities  of  these  initial  cascade  processes  are  im¬ 
mediately  evident.  Ke  formally  integrate  this  expression  to  obtain 


-31- 


where,  in  general,  T  (t)  is  given  by  Eq.  (13).  For  the  specific  case  in 
which  we  regard  I  as  a  constant,  a  change  of  variables  in  Eq.  (18)  is  use¬ 
ful. 


I  I 

When  dt  /dT^  is  substituted  from  Eq.  (9) ,  a  particularly  simple  form  of 
the  electron  density  equation  becomes 


This  integral  may  be  evaluated  in  terms  of  elementary  functions,  and  the 
exponential  integral.  First,  however,  we  define 
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with  these  expressions 


-X 

For  I  >>  kTe  ,  we  have  1/X  <<  1.  Thus,  we  may  neglect  e  °  and  E, 

p  o  o  1 

(X  ) ,  and  we  obtain 
o 


-X 

XQ  A  “x - E1  (X)  (24) 


with  A  independent  of  the  neutral  particle  density,  n  ,  and  inversely 
proportional  to  the  laser  flux,  I.  Thus,  for  times  at  which  the  elec¬ 
tron  temperature  is  notably  higher  than  the  initial  cool  gas,  i.e.. 
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we  have. 


In 


n 


n 


(26) 


For  typical  values  of  fluxes  and  pressures  near  atmospheric,  the  first 
term  will  dominate,  and  the  exponential  term  quickly  approaches  its  asym¬ 
ptotic  value;  thus,  our  final  simple  result  is  characterized  by  a  stronger 
dependence  upon  neutral  density  than  upon  the  laser  flux.  If  we  now  wish 
to  consider  the  detector  responses  of  the  schematic  experiment  in  Fig.  9, 
we  may  make  these  measurements  at  two  locations  simultaneously  since  the 
same  pulse  is  causing  breakdown  at  two  separate  points.  The  maximum  de¬ 
crease  of  the  signal  will  be  a  direct  measurement  of  the  integrated  elec¬ 
tron  density  through  the  plasma,  and  the  time  at  which  the  two  signals  de¬ 
crease  relative  to  one  another  can  be  monitored  accurately. 


(27) 


where  the  subscripts  1  and  2  refer  to  simultaneous  measurements  of  the 
attenuation  of  the  probe  beams  passing  through  the  two  focal  points. 
Thus,  by  monitoring  the  relative  intensities  just  before  opacity  is 
reached,  it  would  seem  that  a  sensitive  measurement  of  the  densities 
(neutral)  at  the  two  different  focal  points  can  be  made.  The  problem 


-34- 


of  fluctuation  of  the  shot  to  shot  laser  pulse  now  poses  no  problem,  and 
the  dependence  upon  the  initial  electron  density  (assuming  threshold  can 
be  reached)  is  weak.  We  should  also  note  that  a  more  complete  analysis 
retaining  loss  terms  may  give  a  somewhat  weaker  density  dependence,  never 
theless,  an  exponential  dependence  is  still  to  be  expected. 

3.1.3  Conclusions 

It  is  well  known  that  a  TEA  laser  pulse  can  ionize  gas  at  some 
distance,  and  that  the  rise  time  of  the  electron  density  is  a  function 
of  several  variables,  the  laser  flux,  the  neutral  density,  the  initial 
electron  density,  and  the  ionization  potential  being  among  the  most  im¬ 
portant.  One  of  the  drawbacks  in  attempting  to  use  such  a  breakdown  to 
obtain  quantitative  information  on  the  state  of  the  gas,  is  that  fluc¬ 
tuations  in  the  power  output  of  the  TEA  laser  may  well  be  larger  than 
the  changes  in  the  state  of  the  gas  from  point  to  point  in  the  flow  field 
With  the  use  of  a  beam  splitter  to  assure  that  at  each  shot,  essentially 
the  same  shape  of  the  pulse  is  causing  the  breakdown,  this  problem  can 
be  avoided.  Our  calculations  indicate  that  this  method  may  also  have  a 
relatively  high  degree  of  sensitivity. 
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3.2  Ultraviolet  -  Visible  Fluorescence  Probes  to  Monitor  Density 

Fluctuations 

There  are  four  potential  sources  of  photoexcitation  which  may  yield 
a  fluorescence  signal  containing  information  about  turbulent  density 
fluctuations.  These  are:  laser-induced  fluorescence,  laser-induced 
photodissociative  excitation,  two-or  three-photon  photodissociative 
excitation,  and  chemiluminescent  reaction  following  photolysis.  In  the 
following,  we  will  examine  each  of  these  options  briefly,  after  setting 
initial  limits  on  reasonable  signals  to  guide  our  choice. 

As  explained  in  Section  4  of  this  report,  to  achieve  10%  precision 
in  the  measurement  of  1%  density  fluctuation,  10^  photons  must  be  detected. 

If  the  detector  system  consists  of  a  collection  lens  at  f/5,  an  inter¬ 
ference  filter  of  50%  transmission  and  a  photomultiplier  of  20%  quantum 

9 

efficiency,  then  4  x  10  photons  must  be  emitted  in  the  region  of  interest  in 
0 

order  to  detect  10  photons.  If  we  neglect  photon  losses  resulting  from 
electronic  quenching,  (an  assumption  which  is  valid  for  all  atmospheric 

quenchers  except  N  ,  O  and  HO  even  at  gas  kinetic  quenching  efficiency), 

^  ^  -13 

and  assume  that  we  are  interested  in  sampling  a  volume  of  2.5  x  10  cm  (a 

volume  roughly  6.3  mm  on  a  side),  then  we  must  have  excited-state  densities 

in  the  viewing  region  of  1.6  x  1010  particles  cm.3  Thus,  any  species  in  the 

atmosphere  whose  abundance  is  less  than  this  can  immediately  be  eliminated 

from  consideration. 

In  order  to  minimize  losses  from  quenching  and  convection  out  of 
the  field  of  view,  candidate  fluorescent  species  should  have  a  radiative 
lifetime  less  than  2ys.  This  requirement  essentially  restricts  the  candidate 
species  to  atoms,  diatomic,  and  perhaps  a  few  triatomic  molecules.  In 
considering  species  in  the  normal  atmosphere,  and  possible  fragments  which 
could  be  produced  from  these  species  in  photolysis,  we  can  find  only  two 
possible  fluorescense  candidates.  Molecular  oxygen  might  be  observed  by 
laser-induced  fluorescence  on  the  Schumann-Runge  bands  in  the  ultraviolet. 
Fluorescence  from  OH  might  be  observable  in  two-photon  photodissociative 
excitation  of  H20. 
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The  output  of  an  ArF  laser  at  193  run  overlaps  the  P  (15,  17,  19) 

and  R  (17,  19,  21)  lines  in  the  v'=4  to  v"=0  transition  of  the  O.B^I  -X  £ 

23  2  u  g 

band.  In  the  thin- target  limit  with  Av  >>Av  ,  ,  the  fractional 

L  absorber 

absorption  of  laser  light  by  any  of  these  lines  is  given  by 

/k  Jtdv  Tre^N&f 

A  =  V  - ^  (2£ 

Av  m  c  Av 

L  e  L 


where  is  the  frequency-dependent  absorption  coefficient,  N  is  the  number 
density  of  the  absorbing  species,  f  is  the  oscillator  strength  of  the  transi¬ 
tion  of  interest,  £  is  the  pathlength  of  the  absorption,  Av  is  the  FWHM  of 

-1  L 
the  laser  line  in  cm  and  the  other  symbols  have  their  usual  meanings.  For 


the  six  lines  of  interest,  we  get 


s  f  band  l 
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6  lines 


where  f  band  is  the  band  oscillator  strength  (3.04  x  10  ^  for  the  (4,  0)  band)^ 

°  18  —3 

N  is  the  number  density  of  oxygen  (4.8  x  10  molecules  cm  at  sea  level), 

2  -1  ' 
i  is  0.1  cm,  Av  is  134  cm  and  the  terms  in  the  sum  reflect  the  Boltzmann 
L 

population  distribution,  and  the  rotational  line  strengths  of  the  lines  under 

consideration.  We  calculate  that  A  =  7.2  x  10."’  If  the  ArF  laser  has  an 
22  12 

output  of  50mJ,  then  3.6  x  10  photons  will  be  absorbed  in  the  viewing 
region. 


Not  all  of  the  absorbed  photons  will  fluoresce,  however.  The 
Schumann-Runge  bands  of  C>2  predissociate,  so  that  only  a  fraction  of  the 
excited  states  produced  in  the  absorption  of  the  laser  light  have  the 
potential  to  fluoresce.  The  fluorescence  efficiency  is  given  by2<^ 
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where  t,.  is 
diss 


the  predissociative  lifetime  given  by 


T 


diss 
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27TcAvj. 

diss 


1.27  x  10 


-12 


sec 


and  T  .is  the  radiative  lifetime  which  is  187  ns  for  v'  =  47  Thus,  in 
rad 

the  absence  of  quenching,  the  detected  fluorescence  signal  will  be  the 
product  of  the  number  of  pnctons  absorbed,  the  fluorescence  efficiency 
against  predissociation,  and  the  system  detection  efficiency.  We  calculate 
that  the  detector  signal  will  be  6  x  105  photons,  giving  a  detection  pre¬ 
cision  of  13%  for  1%  density  fluctuations. 

It  is  to  be  noted  that  the  quenching  of  the  excited  states  will  be 
unimportant  if  the  ArF  laser  is  able  to  saturate  the  transition^  However, 
the  commercially  available  ArF  laser  does  not  produce  saturation.  Therefore, 
the  detection  precision  may  decrease  if  the  quenching  rates  for  the  excited 
states  are  of  the  same  order  as  In  addition,  the  detection  precision 

will  also  decrease  by  approximately  a  factor  of  2  if  the  measurement  is  made 
at  an  altitude  of  10  Km. 

The  excited-state  number  density  produced  from  two-photon 

.  .  ....  ,26 
dissociative  excitation  is  given  by 


(30) 


where  Y  is  the  two-photon  absorption  cross-section,  typically ss 10  cm  W  ,  V 

-2 

is  the  laser  frequency,  I  is  the  laser  intensity  in  W  cm,  is  the  laser 

pulse  length,  and  the  other  symbols  have  their  usual  meanings.  For  an  ArF 

2 

laser  delivering  50mJ  in  20  ns  to  a  1.5mm  target,  we  calculate  that  the 

-4 

excited-state  number  density  produced  will  be  1.5  x  10  of  the  number  density 

of  the  absorbing  species.  Ir  oiuer  to  produce  excited-state  number  densities 
9  -3 

greater  than  10  particles  cm,  the  parent  molecule  must  have  a  number 

12  -3 

density  greater  than  7  x  10  molecules  cm,  or  a  mixing  ratio  greater  than 
—6 

10.  From  Table  1  we  find  that  this  restriction  limits  the  possible  parent 
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molecules  to  N^,  0 ,  H^O,  CO^  and  CH^ .  In  this  list,  the  only  likely 
candidate  for  photodissociative  excitation  would  be  H^O  which  might  produce  OH 
fragments  emitting  around  308  nm. 

In  general,  single  photon  photodissociative  excitation  requires  more 

energy  than  can  be  met  by  sources  which  may  be  propagated  through  air. 

Approximately  4  eV  would  be  required  to  break  a  chemical  bond,  and  another 

3  eV  to  provide  electronic  excitation.  This  amount  of  energy,  7  eV,  can  be 

delivered  only  by  sources  with  wavelengths  shorter  than  175  nm,  i.e.,  beyond 

the  Schumann-Runge  cutoff  in  normal  air.  There  are  some  exceptions  to  this, 

of  course,  such  as  the  photolysis  of  0  at  250  nm  to  yield  electronically 
1  3 

excited  0  (  D) .  We  cannot  think  of  any  exceptions  which  give  fragments  with 
sufficiently  short  radiative  lifetimes  to  be  useful. 

In  the  past  few  years,  much  effort  has  been  expended  searching  for 
efficient  chemiluminescent  reactions  which  produce  visible  photons.  We 
are  unaware  of  any  successes  in  these  endeavors . 

In  summary,  the  only  sources  of  uv-vis  photoexcitation  which 

might  produce  sufficient  light  intensity  to  allow  measurements  in  the 

percent  range  are  laser  induced  fluorescence  of  0^  on  the  Schumann-Runge 

bands,  and  two  photon  photodissociative  excitation  of  H^O  to  produce  OH 

fluorescence.  Both  of  these  techniques  could  be  tried  in  the  laboratory 

28 

with  a  conventional  off-the-shelf  ArF  laser. 
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4 .  TURBULENCE  CONSTRAINTS 

Making  measurement  in  a  turbulent  environment  requires  certain 
spatial  and  temporal  resolution  capability  in  a  diagnostic  technique.  In 
general,  in  a  turbulent  flow,  one  is  interested  in  obtaining  information 
about  the  small  scale  structure  and  the  energy  spectrum.  The  small  scale 
structure  in  a  turbulent  flow  has  a  characteristic  dimension  given  by  the 
Kolmogorov  microscale. 


-i±) 


where  T)  is  the  small  scale  characteristic  dimension,  V  is  the  kinematic 

m 

viscosity  of  the  medium  and  £  is  the  rate  of  turbulent  dissipation  of  energy. 
Since  n  is  proportional  to  the  one  fourth  power  of  £,  the  turbulence  microscale 
for  a  variety  of  flows  of  interest  range  from  0.5  mm  to  1  cm.  For  our  appli¬ 
cation,  the  density  measurements  are  pertinent  to  the  turbulent  flow  around 
a  turret  mounted  on  an  airplane.  It  is  specified  that  the  spatial  resolu¬ 
tion  required  is  approximately  1  cm. 

In  Section  2,  we  have  seen  that  the  most  promising  three-level  gain 
measurement  technique  yields  a  gain  of  0.017  cm  1.  Now  the  gain  is  directly 
proportional  to  the  density  of  the  medium.  Consequently,  to  measure  1% 

-4 

density  fluctuations  we  require  a  gain  measurement  capability  of  1.7  x  10 
cm  1.  In  order  to  see  how  accurately  this  measurement  can  be  made,  we  calculate 
the  signal  to  noise  ratio  in  the  next  subsection.  We  also  note  that  wind 
tunnel  testing  with  air  seeded  with  CO^  will  significantly  increase  the  sensi¬ 
tivity. 

4.1  Signal  to  Noise  Ratio 


In  Section  2  we  have  seen  that  the  energy  constraint  on  the  pump 
laser  requires  that  the  single  pulse  experiment  be  completed  in  a  time 
less  than  a  microsecond.  In  such  a  short  time,  the  noise  contribution 
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due  to  intrinsic  design  of  the  detector  is  small.  The  background  noise 
due  to  300  K  black  body  radiation  does  not  enter  in  the  calculation  of  signal 
to  noise  ratio  because  we  measure  the  change  in  the  signal  and  not  its 
absolute  value. 

Another  source  of  noise  which  can  significantly  reduce  the  magnitude 

of  the  signal  is  the  probe  laser  absorption  in  the  atmosphere.  As  was 

calculated  in  Section  2,  the  absorption  coefficient  at  the  probe  laser 

-4  -1 

frequency  is  7.5  x  10  cm  .  Since  the  path  length  for  the  probe  beam  is 

expected  to  be  approximately  1  meter,  the  probe  laser  absorption  in  the 

-2 

atmosphere  will  be  7.5  x  10  .  Even  5%  fluctuation  in  this  absorption 

from  shot  to  shot  can  conceivably  wipe  out  the  gain  that  we  are  trying  to 
measure . 

This  problem  of  the  atmospheric  absorption  of  the  probe  beam  can  be 
resolved  if  we  trigger  the  scope  connected  to  the  detector  a  few  microseconds 
before  we  fire  the  pump  laser.  This  arrangement  will  establish  a  reference 
level  and  when  the  pump  laser  pulse  arrives,  we  can  just  measure  the  change 
from  this  reference  line.  Since  the  single  pulse  experiment  is  completed  in 
time  less  than  a  microsecond  and  the  characteristic  time  scale  of  the 
turbulence  is  tens  of  microseconds,  there  will  be  no  change  in  this  reference 
line  in  a  few  microseconds.  From  pulse  to  pulse,  however,  this  reference 
line  will  shift,  but  we  can  still  measure  the  gain. 

Fluctuations  in  the  detected  intensity  due  to  photon  statistics  is 
another  source  of  noise.  This  is  the  most  important  noise  source  which 
limits  the  spatial  resolution  of  the  density  measurement  scheme.  The  CW 
diode  laser  to  be  used  to  probe  the  saturated  medium  delivers  0.5  mW  at 
4.32ym.  In  100  nsecs,  which  is  expected  to  be  the  time  to  complete  a  single 

Q 

pulse  experiment,  the  laser  will  deliver  1.1  x  10  photons  to  the  detector. 

9  ^ 

The  statistical  fluctuations  in  detected  signal  is  (1.1  x  10  )  and,  there¬ 
fore,  the  ratio  of  detected  intensity  fluctuations  to  the  initial  intensity 
called  the  noise  level,  is  3.0  x  10  .  This  calculated  noise  level  depends 

upon  the  power  of  the  probe  laser,  the  pulse  duration  of  the  pump  laser  and 
the  effective  vibration  relaxation  rate  of  the  excited  level  in  C02 . 
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In  Fig.  10,  we  show  the  S/N  ratio  as  a  function  of  spatial  resolution 
for  1%  density  fluctuations.  It  can  be  seen  that  in  the  spatial  range  of 
interest  to  us,  1  cm,  the  S/N  ratio  is  fairly  reasonable.  For  instance, 
for  a  path  length  of  1  cm,  it  is  anticipated  that  we  can  measure  1%  density 
fluctuations  to  an  accuracy  of  16%.  It  is  emphasized  here  that  the  S/N  ratio 
as  calculated  above  is  only  approximate  and  at  sea  level.  At  an  altitude  of 
10  Km,  this  ratio  will  drop  by  a  factor  of  3.  However,  if  the  S/N  ratio  is 
not  adequate  to  give  a  reasonably  accurate  measurement,  we  can  always  use  a 
multipass  cell  and  increase  the  effective  path  length  for  gain.  In  Fig.  10 
therefore,  we  have  also  presented  the  S/N  ratio  for  multiple  passes  through 
the  same  saturated  volume.  As  can  be  seen,  a  multipass  cell  allowing  the 
probe  beam  to  pass  several  times  through  the  excited  volume  should  be  sufficient 
to  provide  a  reasonably  accurate  measurement  of  the  density  fluctuations. 

Now  we  calculate  the  temporal  resolution  required  to  make  measurements 
in  a  turbulent  environment  and  see  if  the  three-level  scheme  under  considera¬ 
tion  can  provide  such  temporal  resolution. 

4.2  Temporal  Resolution 

The  highest  frequency  information  which  can  be  obtained  from  a  tur¬ 
bulent  flow  field  is  characterized  by  the  time  scale  of  the  small  scale 
structure.  This  time  scale  is  given  by  the  Kolmogorov  time, 


For  our  application,  the  value  of  the  turbulent  dissipation  rate,  e,  is  not 
known,  but  it  is  believed  that  the  Kolmogorov  time  is  of  the  order  of  tens 
of  microseconds .  This  time  scale  is  longer  than  the  pulse  length  of  the 
pump  laser,  which  is  approximately  six  tenths  of  a  microsecond  as  noted  in 
Section  2.  Thus,  each  sample  is  taken  in  a  time  much  shorter  than  that 
characterizing  the  fastest  flow  motion.  This  means  that  a  particular  sample 
contains  all  the  information  about  the  fastest  flow  motion  at  some  instant  of 
time.  However,  that  is  not  enough, and  to  obtain  frequency  spectrum  informa¬ 
tion  we  need  to  take  samples  at  the  rate  of  tens  of  thousand  samples  per 


SSES 


second.  This  is  because  using  pulsed  devices,  we  can  only  recover  frequency 
information  up  to  the  Hyquist  frequency,  which  is  half  the  frequency  at 
which  the  device  is  pulsed. 

The  commercially  available  high  power  HF  lasers  delivering  power 
levels  of  interest  to  us  can  provide  rep  rates  of  only  1-3  pps .  These  rep 
rates  are  insufficient  to  provide  frequency  spectrum  information.  However, 
if  we  assume  that  the  turbulence  flow  field  is  random  and  isotropic,  i.e., 
taking  samples  at  different  times  is  equivalent  to  taking  samples  at  different 
points  in  space,  we  can  obtain  information  about  the  fastest  flow  motion 
using  the  commercially  available  lasers.  To  see  how  this  is  possible, 
consider  Fig.  11.  The  probe  beam  is  made  to  pass  through  a  beam  splitter 
which  divides  the  beam  into  two  parts.  These  parts  may  not  be  of  equal 
intensity  because  we  are  only  interested  in  the  change  in  intensity.  These 
two  parts  then  serve  as  two  probe  beams  and  sample  the  saturated  gas  at 
two  different  points  in  space.  By  varying  the  distance  between  these  two 
beams,  we  can  determine  the  spatial  autocorrelation  coefficient,  a 

p{x,t)  p(x  +  a) ,  t 

a  =  -  (33) 


where  a  is  the  distance  between  the  two  points  where  the  density  measurement 
is  being  made  at  a  particular  instant  of  time.  If  the  turbulence  is  homo¬ 
geneous,  a  is  only  a  function  of  a,  i.e.  a  =  a  (a) .  By  taking  the  Fourier 
transform  of  a  (a) ,  we  can  obtain  the  energy  spectrum  of  the  fastest  flow 
motion  in  terms  of  wave  numbers . 

From  Eq.  (33) ,  the  power  spectrum  in  terms  of  frequency  can  also 
be  obtained  if  we  assume  that  the  turbulence  is  random.  Then  the  space 
and  time  coordinates  in  Eq.  (33)  can  be  interchanged  and  by  taking  Fourier 
transform  in  frequency  space,  we  can  obtain  the  power  spectrum. 

If  the  turbulence  in  a  flow  field  is  not  completely  random,  we  may 
not  be  able  to  recover  frequency  spectrum  information  but  we  can  always 
generate  probability  distribution  functions.  If  N  is  the  total  number  of 
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density  measurements  made,  and  n  is  the  number  of  measurements  that  lie 
between  p(t)  and  p(t  +  At),  then  the  probability  density,  B  (p)  is  defined 
as 


B  (p)  Ap  =  Lim 

N-t-co 


n_ 

N 


The  probability  density,  B(p)  satisfies  the  following  normalization  condition 


B (p)  >  0  , 


-  OO 


d  p  =  1 
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5.  EXPERIMENTAL  CONSIDERATIONS 


In  the  first  four  sections  of  this  report,  we  have  seen  that  the 
three-level  gain  measurement  technique  which  involves  using  2.76  um  HF 
laser  for  pumping  at  4.32  vim  can  possibly  yield  useful  information  about 
the  density  fluctuations  in  a  turbulent  flow.  We  have  also  seen  that  for 
1%  density  fluctuations,  the  S/N  ratio  is  such  that  a  reasonably  accurate 
measurement  can  be  made  with  spatial  resolution  of  1  cm.  In  this  section, 
we  explore  the  possibility  of  performing  a  simulation  experiment  in  the 
laboratory . 

The  HF  laser  power  in  2  +  1  P  (4)  line  needed  to  saturate  000  +-)'021 

2 

band  in  CC>2  gas  must  be  greater  than  10  MW/cm  ,  while  the  pulse  length  must 
be  less  than  or  equal  to  six  tenths  of  a  microsecond.  Two  types  of  high 
power  pulsed  HF  lasers  are  available.  Both  of  these  models  are  manufactured 
by  Lumonics  Research  Limited.  One  is  the  Model  210,  which  can  deliver  3.6  MW 
of  total  power,  and  the  other  is  the  Model  203,  which  can  deliver  1  MW  of  total 
power.  The  salient  features  of  both  of  these  laser  systems  are  presented  in 
Table  3.  The  Model  210  has  been  discontinued  by  Lumonics  Research  Limited,  but 
can  be  obtained  custom  made  at  an  approximate  cost  of  $40,000.  The  Model  203 
is  a  modified  version  of  Model  103.  At  PSI,  we  have  Model  103  and  it  might  be 
possible  to  modify  this  model  to  obtain  1  MW  of  power  with  HF  gas. 

2 

From  Table  3,  it  can  be  seen  that  greater  than  6  J/cm  fluence  for 

laser  line  2  +  1  P  (4)  can  be  obtained  if  the  beam  is  focussed  to  less  than 

1.0  mm  diameter  for  Model  210  and  0.5  mm  diameter  for  Model  203.  Thus,  the 

2 

available  commercial  lasers  cannot  saturate  a  path  greater  than  1  mm  in 
cross-section.  To  obtain  a  gain  length  of  1  cm,  one  may  proceed  in  two  ways. 
First,  in  the  laboratory,  the  scheme  could  be  tested  at  a  reduced  pressure. 

For  C02  gas,  the  collision  broadened  regime  extends  down  to  10  torr  of  pressure. 
In  a  collision  broadened  regime,  the  saturation  intensity  is  directly  pro¬ 
portional  to  the  pressure  through  the  vibrational  relaxation  rate  constant. 
Consequently,  by  testing  the  scheme  in  a  wind  tunnel  at  a  reduced  pressure, 
say  10  torr  (keeping  the  concentration  of  C02  unchanged^  we  can  relax  the 
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TABLE  3 


Relevant  Information  on  Commercially  Available  High  Power  HF  Lasers 


Pump  Laser 


Type 

Lumonics  -  210A 

Lumonics  -  203 

Power  (multiline) 

3.  6  MW 

1  MW 

Power  (single  line) 

0.  6  MW 

100  KW 

Power  (2  -»  IP  (4)  line) 

0.  15  MW 

25  KW 

2 

Beam  dia  to  obtain  6j/cm 

1.  0  mm 

0.  50  mm 

Pulse  width 

0.  34  psec 

0.  5  psec 

Cost 

$  35  K  -  40  K 

$  22  K 

requirement  on  saturation  fluence  by  a  factor  of  76.  This  arrangement  will 
permit  us  to  obtain  a  gain  length  of  approximately  8.7  mm  with  the  Model  210 
laser. 

The  second  possibility  is  to  use  a  cylindrical  focussing  lens  which 
will  permit  us  to  obtain  rectangular  cross-sections  in  the  focal  plane. 
Focussing  the  beam  to  0.1  mm  in  one  dimension,  we  can  obtain  a  gain  length 
of  1  cm  in  the  other  dimension.  This  seems  a  more  promising  method  for 
obtaining  saturation  and  a  cm  beam  width. 

The  laser  to  be  used  to  measure  gain  at  4.32  ym  is  the  CW  diode 
laser  manufactured  by  Laser  Analytics  Inc.  This  laser  can  be  fine  tuned  to 
provide  a  laser  frequency  which  is  in  resonance  with  the  P(16)  line  of 
021  ■*-*■  020  band  in  CO^.  The  laser  operates  at  cryogenic  temperatures  and 
requires  a  lens  assembly  to  collimate  the  beam.  Definite  information 
about  the  amplitude  stability  of  this  laser  over  a  microsecond  time  scale 
is  not  available,  but  it  is  believed  that  the  amplitude  stability  is  better 
than  one  part  in  one  hundred  thousand  during  the  time  of  the  pumping  pulse. 

If  it  is  necessary  to  use  a  multipass  cell  to  increase  the  effective 
gain  length,  this  can  also  be  obtained  from  Laser  Analytics  Inc.  This  multi¬ 
pass  cell  contains  two  confocal  mirrors  placed  at  two  ends  of  a  one  meter  cell 

By  appropriately  varying  the  focal  lengths  of  the  two  mirrors,  we  can  adjust 

.  .  2 
the  beam  cross-section  (in  the  focal  plane)  anywhere  in  the  range  1  mm  to 

2 

1  cm.  Use  of  a  multipass  cell  requires  transfer  optics  and  it  is  apparent 
that  the  cost  of  testing  the  technique  will  increase. 

One  advantage  of  making  the  gain  measurement  at  4.32  ym  is  that 

sensitive  detectors  are  available  in  this  range  of  the  infrared.  Indium 

antimonide  detectors  can  detect  4.32  ym  radiation.  These  detectors  have  a 

rise  time  of  0.2  ysec  and  a  D*  of  lO1^  cm  H  ^  watt  1.  The  noise  equivalent 
-9  2 

power  is  10  watts,  which  is  small  in  comparison  to  the  power  of  CW  diode 
lasers . 

At  PSI,  the  density  fluctuation  diagnostic  can  be  assembled  and 
experiments  can  be  performed  in  gas  cells  at  different  pressures  to  establish 
gain  levels,  signal  to  noise  ratio,  and  determine  if  the  pump  laser  is 


saturating  the  CO^  gas.  Saturation  of  the  gas  is  insured  if  the  small 
signal  gain  is  directly  proportional  to  the  CO^  number  density.  Then  low 
speed  turbulent  flow  experiments  (grid  turbulence)  should  be  performed. 

A  transonic  wind  tunnel  with  a  6"  by  6"  test  section  is  located  in  the 
Division  of  Engineering  at  Brown  University,  and  small  scale  turbulence 
experiments  could  be  performed  using  this  tunnel.  Once  these  experiments 
are  completed  and  the  features  of  this  measurement  technique  have  been 
validated,  large  scale  wind  tunnel  testing  should  be  performed  at  either 
NASA  or  Air  Force  laboratories. 
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6.  CONCLUSIONS 


In  this  report,  several  three-level  gain  measurement  schemes  have 
been  examined  from  the  point  of  view  of  monitoring  density  fluctuations 
in  a  turbulent  environment  in  the  atmosphere .  Two  other  techniques ,  that 
of  monitoring  laser  induced  fluorescence  signal  and  the  use  of  gas  break¬ 
down  as  a  diagnostic  tool  have  also  been  examined.  Our  analysis  shows 
that  out  of  all  the  schemes  considered  in  this  report,  the  three-level 
scheme  which  involves  using  a  2.7  ym  HF  laser  to  saturate  000  -*-►  021 
band  in  CO ^  and  gain  measurement  at  10.6  ym  with  a  CW  diode  laser  is  the 

most  promising  scheme.  The  pump  laser  fluence  required  for  saturation 

2 

comes  out  to  be  equal  to  or  greater  than  6  J/cm  which  is  achievable  with 
commercially  available  lasers.  Quantitative  information  about  the  energy 
spectrum  of  the  density  fluctuations  can  be  obtained  by  making  simultaneous 
measureme'nts  at  two  different  points  in  space.  Furthermore,  remote  and 
spatially  resolved  measurements  providing  information  about  the  fastest 
flow  motion  can  be  made  with  reasonable  degree  of  accuracy. 

Another  technique  which  has  shown  some  promise  of  being  able  to 
measure  1%  density  fluctuations  in  a  turbulent  environment  is  the  monitoring 
of  laser  induced  fluorescence  of  on  the  Schumann-Runge  bands.  It  is 
recommended  that  both  the  schemes  mentioned  above  should  be  tried  in  the 
laboratory  prior  to  field  experiments. 

We  wish  to  emphasize  that  even  though  the  various  schemes  in  this 
report  are  considered  from  the  point  of  view  of  monitoring  atmospheric 
density  fluctuations,  the  three-level  scheme,  in  principle,  can  also  be 
used  to  measure  species  concentration  in  a  variety  of  other  applications . 

For  instance,  there  is  considerable  interest  in  determining  the  extent  of 
spatial  non-uniformities  in  gain  in  high  power  HF  lasers.  The  extension 
of  the  three-level  scheme  to  this  application  is  straightforward.  Other 
applications  where  the  three-level  scheme  may  be  of  interest  include  the 
measurement  of  species  concentration  in  MHD  flow  channels  and  in  turbulent 
combustion  flows. 
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APPENDIX  I 


CALCULATION  OF  SATURATION  INTENSITY 


We  define  the  saturation  intensity  for  a  given  transition  as  that 
intensity  which  results  in  zero  gain  (in  the  limiting  sense  only)  on  that 
transition.  In  a  two-level  vibration  -  rotation  transition,  the  small  signal 
gain  is  given  by  Eq.  (3).  Zero  gain  implies  that 


g2,J 
91,J  +  a 


0 


or 


’2,J 


2,J 


*1,J  +  a 


1 ,  J  +  a 


(AX-1) 


where  a  =  +1,  0,  -1  corresponds  to  P,  Q  and  R  branch  transitions. 

When  Eq.  (AI-1)  is  satisified,  absorptions  and  stimulated  emission 
balance  one  another  and  no  further  excitation  results  with  increased  pump 
intensity.  Eq.  (AI-1)  also  implies  that  if  level  1  is  the  ground  vibrational 
level,  the  collisional  and  radiative  relaxation  of  excited  level  (2,J)  other 
than  via  stimulated  emission  is  negligible. 


Now  we  derive  an  expression  for  the  laser  intensity  required  to 
saturate  a  given  rotation  -  vibration  transition.  We  consider  vibrational 
levels  1  and  2.  Within  each  vibrational  manifold  is  a  large  number  of 
rotational  levels.  The  rate  equations  for  levels  (2,J)  and  (1,  J  +  a)  may 
be  written  as 


-Kn  -  R  (n.  ,  -  n.  _  'j  + 

2, J  V  2 , J  1 , J  +  a / 
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(AI-3) 


Here  K  is  the  effective  rate  constant  for  vibrational  relaxation  of  level 

(2,J)  via  processes  other  than  stimulated  emission,  R  is  the  stimulated 

emission  rate  constant  (and  it  is  assumed  that  the  laser  is  resonant  with 

only  one  rotational  line),  and  a  is  the  rotational  relaxation  rate  constant. 

r 

For  simplicity,  in  Eqs .  (AI-2)  and  (AI-3)  we  have  taken  degeneracy  ratio 
equal  to  1. 

The  stimulated  emission  rate  constant,  R,  is  proportional  to  the 
Einstein  B  coefficient: 


R  = 


1,J 


+  a;2,J 


(AI-4) 


where  g  (v,  Vq)  is  the  line  shape  function  and  is  the  line  center  frequency. 

Summation  of  Eqs.  (AI-2)  and  (AI-3)  over  all  rotational  states  yields 
the  rate  equations  for  vibration  level  2  and  1. 


(AI-5) 


■  R  ("j  -  Bl)  *  K  «2 


(AI-6) 


Analytical  solutions  of  Eqs.  (AI-2)  through  (AI-6)  for  population 


inversion  cannot  be  obtained.  However,  if  we  assume  that  R,  »  K 
we  can  solve  for  (N,  -  N  ) . 


=  Boltzmann  factor  = 


The  slowest  rate  for  population  build  up  in  vibrational  level  2,  a2  can  be 
simplified  as  follows: 
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Thus,  in  the  absence  of  any  vibrational  relaxation,  the  steady  state  in 
population  will  be  obtained  when  2BfR  «0^.  Hence,  the  condition  of 
population  saturation  which  yield  a  zero  gain  on  vibration  levels  1  and  2 
and  on  any  of  the  rotational  lines  associated  with  these  two  vibration 
levels  is 


a2  »  2BfR  »  K  (AI-9) 

Comparing  this  condition  with  the  saturation  condition  for  a  two  level  atomic 
system  (2R  »  K) ,  we  note  that  if  the  rotational  relaxation  is  allowed  to 
occur  during  the  time  when  the  laser  radiation  is  pumping  molecules,  the 
stimulated  emission  rate  must  be  increased  by  a  Boltzmann  factor.  Thus  the 
intensity  requirement  to  saturate  two  vibrational  levels  when  only  one 
rotation  line  is  resonant  with  the  laser  radiation  is  much  more  stringent. 

It  is  to  be  noted  that  when  2  B^R  <  K,  one  can  still  obtain  a  steady 
state  population  distribution,  but  in  this  case,  the  gain  between  levels  2 
and  1  will  not  be  zero. 

Now  we  check  the  possibility  of  satisfying  Eq.  (AI-9)  for  molecules 

6  —I 

of  interest  to  us.  For  CO  at  1  atmosphere,  K  w  1.6  x  10  sec  while 
9-1  2 

0^  —  6.6  x  10  sec  .  Thus  K  and  0^  differ  by  three  orders  of  magnitude 
and  one  can  conveniently  choose  a  suitable  value  of  the  laser  intensity  which 
will  satisfy  Eq.  (AI-9) . 

We  can  calculate  the  lower  bound  on  the  laser  intensity  needed  to 
saturate  a  vibration  -  rotation  transition  (Ig)  from  the  following  equation: 


2  R  Bf  =  K  . 


Making  use  of  Eq.  (AI-4)  and  assuming  that  the  laser  radiation  is  a  delta 
function  in  frequency,  we  get 


2  B  B 

f  1 ,  J  +  a ;  2 ,  J 


K 
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For  a  pressure  broadened  transition,  the  line  shape  function  is 
Thus , 


I 

s 


16  TT  2  Kh  VQ3 
2  Bf  c  A 


Lorentzian . 


(AI-10) 


where  A  is  the  Einstein  coefficient  for  spontaneous  emission. 


APPENDIX  II 


INCREASED  INVERSION  EFFICIENCY  IN  OPTICALLY* 
PUMPED  TWO -LEVEL  MOLECULAR  SYSTEMS 


ABSTRACT 

It  is  shown  how,  under  certain  circumstances,  optically  pumping  an 
R  branch  transition  at  high  J  number  in  a  molecular  gas  may  result  in  a 
situation  in  which  over  half  of  the  particles  are  in  the  excited  vibrational 
manifold.  This  can  lead  to  increased  gain  in  transitions  to  intermediate 
lying  states  and  increased  fluorescence  signal  in  laser  induced  fluorescence 
experiments.  In  addition,  gain  can  also  be  obtained  on  pumping  band  in 
transitions  which  have  frequencies  smaller  than  the  pump  laser  frequency. 
Such  an  arrangement  can  be  employed  to  obtain  a  wavelength  shift. 


*To  be  published  in  Optics  Letters,  December  1979. 
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The  purpose  of  this  note  is  to  specify  conditions  in  a  molecular  gas 

under  which  more  than  half  the  particles  in  the  ground  vibrational  state  may 

be  optically  excited  to  a  higher  vibrational  manifold.  In  a  variety  of 

applications  there  is  strong  interest  in  maximizing  the  excited  state 

29 

number  density:  nonlinear  saturation  spectroscopy  with  lasers,  generation 

of  laser  oscillations  on  coupled  transitions  in  an  optically  pumped  three- 

level  system30  33  and  monitoring  gas  density  using  laser  induced  fluorescence 
6 

in  a  gas  flow.  In  a  two-level  atomic  system,  optical  pumping  can  result  in 
the  excitation  of  at  most  half  the  lower  energy  state  particles  (assuming 
that  the  degeneracy  of  the  two  levels  coupled  to  the  laser  radiation  is  equal) . 
Once  the  number  of  particles  in  the  lower  and  upper  level  becomes  equal, 
saturation  sets  in  and  no  further  excitation  results  from  increased  pump 
intensity.  In  a  two-level  molecular  system,  however,  as  a  consequence  of 
competing  relaxation  processes,  it  will  be  shown  how  more  than  half  of  the 
lower  vibrational  state  particles  may  be  excited  to  a  higher  vibrational 
level . 

We  consider  a  two-level  molecular  system  with  vibrational  levels  1 

and  2.  Within  each  vibrational  envelope  is  a  large  number  of  rotational 

levels.  We  assume,  for  simplicity,  that  the  pumping  laser  frequency  and 

line  width  are  such  that  absorption  occurs  only  on  one  rotational  line.  In 

a  molecular  system,  typically,  the  rotational  relaxation  rate  is  2  -  3  orders 

of  magnitude  larger  than  V-V  or  V-T  energy  transfer  rates.  Consequently, 

we  can  obtain,  by  a  proper  choice  of  laser  intensity,  a  characteristic 

pumping  time  (T  )  which  is  much  larger  than  T  but  much  less  than  T  ,  where 
p  R  V 

Tr,  and  Ty,  are  the  characteristic  times  for  rotational  and  vibrational 
relaxation,  respectively.  Under  these  conditions,  if  the  pumping  laser  is 
tuned  to  one  rotational  line,  and  if  saturation  occurs,  the  population  of 
these  two  upper  and  lower  vibrational-rotational  levels  will  be  equal . 

However,  this  condition  does  not  necessarily  imply  that  only  half  of  the 
total  molecules  in  both  vibrational  states  can  be  excited  to  the  higher 
vibrational  state. 
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We  consider  the  following.  Let  the  upper  vibrational-rotational 
level  coupled  to  the  radiation  be  denoted  by  (2,J) .  Then  the  lower  vibra¬ 
tion-rotational  level  is  (l,J+a) ,  where  a  =  +1,  0,  and  -1,  respectively, 
corresponds  to  a  P,  Q  or  R  branch  transition.  In  the  presence  of  a  strong 
saturating  field,  the  number  of  molecules  in  level  (2,J)  is  given  by 


n 


2,  J 


g 

g 


2 ,  J 
1,  J+a 


‘1, 


J+a 


(AII-1) 


where  g  j  and  g1>J+a  are  the  degeneracies  of  levels  (2 , J)  and  (1,  J+a) 

respectively.  Since  XR  «  T^,  complete  rotational  equilibrium  exists  in  the 

vibrational  levels  1  and  2.  Consequently,  the  number  of  molecules  n  and 

2 ,  J 

nl  JHa  are  9^-ven  by  the  Boltzmann  distribution 


hcB2  r 

n2,J  =  N2  ~kT~  (2J  +  exp  f  B2hcJ(J  +  DAT 


] 


(AII-2) 


hcBi  r  i  r 

nl, J+a  =  N1  “kT-  2  (J+a)  +  1  exp  r  B1hc(J+2)  ( J+a+1) /kT 


(AII-3) 


where  N^,  N2  are  the  total  number  of  molecules  in  the  vibrational  levels 
1  and  2  and  Bi  is  the  rotational  constant.  From  Eqs .  (AII-1),  (AII-2)  and 
(AII-3),  the  ratio  of  total  population  in  vibrational  manifolds  1  and  2, 
Nj/N^  can  be  determined  as  follows : 


N 

2  _  exp  (-  a(J  +  a) ) 

N1  exp  (-  a(J) ) 


(AII-4) 


where 


a(J) 


BhcJ  (J+l) 
kT 


(AII-5) 


in  Eqs.  (AII-4)  and  (AII-5),  for  simplicity  we  have  assumed  that  the  rotational 
constants  for  levels  1  and  2  are  equal,  i.e.,  =  B^  =  B. 

If  the  total  number  of  molecules  is  N  (N  =  N,  +  N  ) ,  we  can  calculate 

o  o  1  2 

VNo  from  Eq.  (AII-4)  as 

^2  _  exp  (-ot(J+a) _ 

exp  (-a  ( J+a) )  +  exp  (-ct(J))  ^ 

In  Figure  12,  the  ratio  N  /N  as  function  of  the  Boltzmann 

factor,  BF  (B^  =  j/N^)  ^or  ^ '  2  and  R  branch  pumping  in  HF  .and  CO.  These 

are  merely  representative  cases  and  the  conclusions  drawn  in  the  following 

pertains  to  any  optically  pumped  two* level  molecular  system.  It  can  be  seen 

that  for  both  hF  and  CO  molecules,  more  than  half  the  particles  can  be  excited 

to  vibrational  level  2  if  the  pump  laser  is  tuned  to  a  high  J  rotational  line 

on  a.i  R-brarch  transition.  In  fact,  if  we  pump  far  enough  in  the  tail  of  the 

rotational  distribution  on  an  R-branch  transition,  we  can  virtually  empty  the 

ground  state.  The  important  parameter  which  determines  the  extent  of  excita- 

_ 

tion  to  level  2  is  the  rotational  partition  function,  Q  (Q  =  r— — )  .  If  Q  is 

r  r  hcB  *r 

small,  for  the  same  BF,  one  can  excite  more  particles  to  the  excited  level  2 
than  when  Q^_  is  large.  Since  is  directly  proportional  to  the  temperature. 
Fig.  12  also  shows  that  it  is  advantageous  to  pump  a  gas  at  lower  temperature 
in  order  to  increase  the  excited  state  number  density. 

To  consider  gain  on  a  transition  between  level  (2,J)  and  a  level 
intermediate  between  1  and  2  (and  essentially  unpopulated  thermally) ,  we 
need  only  multiply,  Eq.  (AII-6)  by  a  Boltzmann  distribution.  We  also  examine 
the  possibility  of  gain  between  vibrational  levels  1  and  2.  The  population 
inversion 
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at  300K 


actor 


AN, 


*2,J‘ 


2,J 


'  gl,  J'+a'  1'J'+a'  /Nc 


where  a'  =  +1,  o,  -1,  for  P,  Q  and  R-branch  transitions  respectively,  is 
given  as  follows: 


AN  r  n 

J  _  (2  J'-t-l)  Bhc  exp  (-g(J))  exp  (-g(J'+a'))  x 

Nq  kT  [exp  (-a(J+a))  +  exp  (-a(J))] 


“ 

- 


exp 


|^2aJ  -  2a' J'  +  a  -  a' 


+  a 


(AII-7) 


In  a  thermalized  rotational  distribution,  the  population  inversion  on  a  P- 

branch  transition  is  larger  than  on  Q-  or  R-branch  transition.  Consequently, 

for  a  given  rotational  partition  function,  the  maximum  population  inversion 

is  obtained  when  the  rotational  level  with  the  most  particles  in  the  excited 

states  (J1  )  is  involved  in  P-branch  transition.  In  Fig.  13  we  show  this 

max 

maximum  population  inversion  normalized  by  the  total  number  density  as  function 
of  for  P,  Q,  and  R-branch  punping.  It  can  be  seen  that  a  very  large  popula¬ 
tion  inversion  is  predicted  when  pumping  is  done  or.  R-branch  transition.  One 
interesting  result  to  be  noted  from  Fig.  13  is  that  even  though  P-branch  pumping 
excites  less  than  half  the  ground  state  particles  to  the  vibration  level  2,  one 
can  still  obtain  population  inversions  on  P-branch  transitions  having  frequency 
smaller  than  that  of  the  pump  laser.  This  population  inversion,  however,  is 
small  and  is  due  to  the  fact  that  the  frequency  of  P-branch  transition  decreases 
as  J  is  increased. 

It  is  to  be  noted  that  if  the  pumping  laser  is  tuned  farther  and  farther 
away  from  the  maxima  of  the  rotational  distribution,  the  number  of  particles 
available  in  one  rotational  level  keeps  on  decreasing.  Consequently,  to 


-62- 


Fig.  13  Normalized  maximum  population  as  function  of  B 


cause  saturation,  the  pumping  intensity  must  be  increased  (Saturation 
Intensity  ~  absorption  coefficient  S .  However,  in  the  collision 
broadened  regime,  the  saturation  intensity  is  directly  proportional  to 
the  pressure  through  the  characteristic  time  for  vibrational  relaxation. 
Therefore,  in  the  collision  broadened  regime,  by  reducing  pressure  it 
is  possible  to  cause  saturation  at  reduced  pumping  intensity. 

The  extreme  depletion  of  the  ground  state  as  predicted  above  may 
be  pertinent  to  studies  in  the  reaction  kinetics  of  excited  states. 
Increased  population  in  the  excited  state  can  also  lead  to  an  increase  in 
the  fluorescence  signal  in  laser  induced  fluorescence  experiments'*  and 
increased  gain  in  transitions  to  intermediate  lying  states.  Furthermore, 
pumping  on  high  J  values  in  R-branch  transitions  and  obtaining  lasing  on 
other  transitions  which  have  inverted  populations  can  provide  us  with  an 
efficient  method  of  frequency  shifting.  This  can  be  applied  to  the  selec¬ 
tion  of  lines  that  propagate  through  the  atmosphere  with  less  attenuation, 
and  can  also  be  used  to  shift  lines  to  obtain  greater  absorption  in  optical 
pumping  techniques. 
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